Timing analysis of two-electron photoemission 
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We predict a significant delay of two-electron photoemission from the helium atom after absorp- 
tion of an attosecond XUV pulse. We establish this delay by solving the time dependent Schrodinger 
equation and by subsequent tracing the field-free evolution of the two-electron wave packet. This 
delay can also be related to the energy derivative of the phase of the complex double photoioniza- 
tion (DPI) amplitude which we evaluate by the convergent close-coupling method. Our observations 
prompt future attosecond streaking experiments on DPI of He which can elucidate various mecha- 
nisms of this strongly correlated ionization process. 
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The attosecond streaking has made experimentally ac- 
cessible the characteristic timescale of electron motion 
in atoms [1, 2]. Recent applications of this technique to 
atomic photoionization, both in the near-infrared (NIR) 
[3] and extreme ultraviolet (XUV) [4] spectral energy 
range, revealed a noticeable time delay between subject- 
ing an atom to a short laser pulse and subsequent emis- 
sion of the photoelectron. While in the NIR photon en- 
ergy range such a delay can be related to nonadiabatic 
tunneling [5], the XUV delay can be, at least partially, 
attributed to the energy dependent phase of the com- 
plex photoionization amplitude [6, 7]. This observation 
is particularly important as it allows for a complete char- 
acterization of the photoionization process in a so-called 
complete photoionization experiment [8] . 

In the case of double photoionization (DPI), all the 
essential information on the many-electron dynamics of 
this strongly correlated ionization process is contained in 
a pair of symmetrized gerade and ungerade amplitudes. 
The moduli of these amplitudes and their relative phase 
can now be determined experimentally [9, 10]. In this 
Letter, we demonstrate that an additional information 
on the individual phases of the DPI amplitudes can be 
supplemented by an XUV time delay measurement. In 
our demonstration, we consider the helium atom driven 
by an XUV attosecond pulse with the same parameters 
as employed in the attosecond streaking experiment on 
Ne by Schultze et al [4]. By solving the time dependent 
Schrodinger equation (TDSE) and by tracing subsequent 
field- free evolution of the two-electron wave packet, we 
establish the apparent "time zero" when each of the two 
photoelectrons leaves the atom. This time depends sen- 
sitively on the photon energy and the energy sharing be- 
tween the photoelectrons. 

To facilitate an individual attosecond streaking in 
a two-electron ionization process, it was suggested in 
Ref [11] to direct emitted electrons parallel and perpen- 
dicular to the NIR field to respectively maximize and 



minimize its streaking effect. We adopt this configura- 
tion and direct two photoelectrons perpendicular to each 
other fci _L fc2- We also distinguish the reference photo- 
electron, the one which would be streaked, and its spec- 
tator counterpart, which influences the reference photo- 
electron via their mutual Coulomb interaction. For sim- 
plicity of our analysis, both photoelectrons are kept in 
the same xz plane with the polarization vector of the 
XUV radiation directed along the z axis. 

The time-dependent calculation of DPI of He was per- 
formed by radial grid integration of the TDSE using the 
Arnoldi-Lanczos method [12]. We used the linearly po- 
larized XUV pulse £{t) — Eq g{t) coswi with the envelope 
function g{t) — cos^(7rf/2Ti) centered at t = 0, which we 
take as the physical "time zero" . The peak field strength 
was Eq = 0.1 a.u., carrier frequency uj = 106 eV and 
T = 27r/w = 39 as. The pulse was turned off outside the 
interval ±Ti, where Ti = AT. 

The field-free solution of the TDSE for t > Ti was 
used to extract information about the DPI process with 
asymptotic photoelectron momenta fci, fc2. This task 
was achieved by tracing time evolution of the wave packet 
state ^i{ri,r2,t) = Pfcj_fe2^'(rJ, rj, i), where the kernel 
of the projection operator was constructed as 



yi,r2\Pki.M\ri,r2) 



*Q.(n)*q.(^2) (1) 



Here '^^.{ri) are one-electron scattering states with the 
ingoing boundary condition describing a photoelectron 
moving in the Coulomb field with Z — 2. The integra- 
tion region is defined as Q = Cli^Q2, where Qi and 
fl2 are spheres in momentum space centered around the 
momentum vectors fei,fe2 so that \qi — ki\ < 0.25fci. 

The wavepacket state ^'i(i) can be expanded over the 
set of the double continua states of the He atom as 



, (2) 
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'^i[ri,r2,t)^jdqidq2f[qi,q2)'iq^,q^{ri,r2)e' 

where E — qi/2 + (7I/2. When both ri,r2 and ri2 are 
large, 

*qi,q2('^i''^2) oc exp[i(qi • n +92 • T-2 +7)] , (3) 



where 7 is the Redmond logarithmic phase [13]. This 
leads to the following asymptotic expression 

^i{ri,2^ oo,t>Ti) ^ dqidq2\f{qi,q2)\ (4) 
X exp|i[arg/(qi,q2) + Qiri + 92^2 + 7 - -Et]| 

The center of the wave packet moves in such a way that 
its phase is stationary with respect to both qi and q2 at 
the points fci,fc2 of the center of the wavepacket: 

(rf/dfe,)[arg/(fei,fc2) + fciri + fc2r2+7--E<] = , 

. (5) 

where ki is either fci or fc2- This gives asymptotic equa- 
tions for the electron trajectories: 



ri x fcj [t - d'j/dEi - daigf{ki,k2)/dEi\ 



(6) 



The term containing derivative of 7 gives logarithmic 
(with t) corrections to the electron trajectory, which can 
thus be represented asymptotically as 



ri{t) - ht - rii{t) X kit^i 



(7) 



Here ioi = dargf{ki,k2)/dEi are the time delays and 
rii(t) are known functions which vary logarithmically 
with t. 

We apply these asymptotic formulas to describe time 
evolution of the maxima of the electron density defined 
asp{r,t) = /|*i(ri,2,i)P((5(r-ri) + 5(r-r2j) dridr2 
and ^'i(ri_2, t) provided by the TDSE solution. As an il- 
lustration of our technique, we consider a DPI process in 
which one photoelectron escapes with energy 8 eV along 
the z-axis and the other with energy 20 eV along the 
X-axis. The snapshot of the electron density in the xz 
plane corresponding to the moment of time t = lAT (10 
field cycles elapsed after the end of the XUV pulse) is 
shown on the top panel of Fig. 1. The figure exhibits 
clearly the two well formed maxima corresponding to the 
center of the wavepacket '^i{t) propagating in the x and 
z directions. A sequence of such snapshots is taken with 
an interval of 2T and the maxima of the electron density 
are traced in time. This procedure defines the trajecto- 
ries ri{t) for both photoelectrons which are exhibited in 
the middle panel of Fig. 2. The free propagation is vi- 
sualized to the straight lines V^t and Vxt. Knowing the 
trajectories, we can compute the left-hand side of Eq. (7) 
which is plotted on the bottom panel of the Fig. 1. At 
sufficiently large time, this quantity should approach the 
constant values of the time delays toi for both photoelec- 
trons. These values are 107 as and 28 as for the 8 eV 
and 20 eV photoelectrons, respectively. We ran an anal- 
ogous simulation with the directions of the fast and slow 
photoelectrons being swapped and obtained very similar 
time delay figures. We also explored the case of the equal 
energy sharing of both photoelectrons and obtained the 
time delay of about 55 as for the 14 eV photoelectrons. 

To relate the time delay to the phases of the DPI am- 
plitudes, we employ the lowest order, with respect to the 
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FIG. 1: Time evolution of the two-electron wave packet for 
8 eV and 20 eV photoelectrons propagating along the z and 
X axes, respectively. Top: the electron density plot in the xz 
plane ai t — lAT. Middle: trajectories of both photoelectrons 
as functions of time measured in numbers of field periods. 
The straight lines visualize the free propagating V^t and V^t. 
Bottom: the effective time delay computed from the LHS of 
Eq. (7). 

field, perturbation theory (LOPT) on the basis of chan- 
nel states of the convergent close-coupling (CCC) method 
[14]. This basis is composed of the products of a Coulomb 
wave Xk (fo'' faster of the two photoelectrons) and a 
positive energy pseudostate (for the slower photoelec- 
tron) which is obtained by diagonalizing the He"*" Hamil- 
tonian in a truncated Laguerre basis. Thus, we write 

vl'(ri,r2,i) {Xk<l>^\mo) Xfe(ri) Mr2) 



(8) 



Here £{u>) = e*"'^£(r) dr is the Fourier transform of 
the XUV field, D is the two-electron dipole operator, 



Eik = + Ci and > > 0. By projecting 

the positive energy pseudostate onto the matching energy 
Coulomb wave fc|/2 = e„2i2: we restore the continuum 
normahzation and phase and write a two-electron wave 
packet in the form of a partial wave expansion: 



*i(n,T-2,0 



dki Di^i^{ki,k2) 



(9) 



y.Rk,iM)Ru,i,{r2)yr'{fi,f2)e-'^'£{E ^ Eo) 

Here 3^/^'^ is a bipolar harmonic and i? — fcJ/2-|-fc|/2. 
The two-electron dipole matrix element is defined as 

>^ Di,i,{kin2){l2k2\\l2n2) (10) 

with Di-^i^{kin2) to be found by integrating the bare 
dipole matrix element with the half on-shell T-matrix 
[14]. 

Knowing the asyniptotics of the radial orbitals Rki oc 
sin [fcr -I- ^i(fc) -f l/fcln(2fcr) — /7r/2] and applying the 
usual saddle-point approximation, we arrive to Eq. (5) 
with the following definition of the DPI amplitude: 



/(fci, fc2) = Di,iAki,k2)yt'' (fci, fc2) 



(11) 



To obtain the time delay toi = dsLTgf{ki,k2)/dEi, we 
have to take the derivative over the energy of the corre- 
sponding photoelectron. However, these derivatives are 
only defined for Ei > E/2 and E2 < E/2. So neither 
of photoelectrons can serve as the reference electron in 
the full excess energy range. In the following, we choose 
the photoelectron 2 to be the reference one. We compute 
the time delay to2 which is defined for E2 < E/2, and 
continue it analytically past the mid excess energy point 
El > E/2. 

An example of the phase plot of the DPI amplitude 
is given in Fig. 2. In the xz plane, geometry of the 
two-electron escape is fully defined by the two azimuthal 
angles 6'i,02- On the top panel of Fig. 2, the phase 
aTgf{Ei,E2) is plotted for the slow electron energies 
i?2 = 2, 4, 6, 10 and 15 eV and variable angle 62 whereas 
the fast electron energy Ei = 40 eV and its angle 61 — 0° 
are fixed. On the bottom panel, directions of the pho- 
toelectrons are swapped and the angle of the slow pho- 
toelectron is fixed at 62 — 0° while the angle of the fast 
photoelectron 9i varies. The phases of the DPI ampli- 
tudes displayed in Fig. 2 can be used to obtain the tim- 
ing information of the two complementary processes in 
which the slow reference photoelectron is directed along 
with (bottom panel) and perpendicular to (top panel) 
the XUV field. 

From inspection of Fig. 2, we see that the phases of 
the DPI amplitudes depend sensitively on the mutual 
photoelectron orientation. However, the spacing between 
the various E2 phase curves does not change significantly 
with the variable photoelectron angle ^2- So the energy 
derivative of the phase of the DPI amplitude, and hence 
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FIG. 2: The phases of the DPI amplitude f{Ei,E2) are plot- 
ted for E2 = 2, 4, 6, 10, 15 eV and Ei = 40 eV. On the top 
panel, the fast photoelectron angle 9i — 0° is fixed and the 
slow electron angle 62 is variable. On the bottom panel it is 
vise versa. 

the effective time delay, does not change very much with 
the relative orientation of the photoelectrons. 
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FIG. 3: Top: cumulative phase plot arg £2) for various 
combinations of the streaked and spectator electron energies 
taken at the perpendicular orientation 612 = 90°. Bottom: 
The time delay to2 = d'AXg f{k\,k2) / dE2 for various parallel 
and perpendicular electron energies. The values of time delay 
obtained from the solution of TDSE are marked by the points 
colored respectively to match the combination of Ei and E2 
from the CCC calculation. 

The vertical line on the both panels of Fig. 2 marks 
the mutual angle of the photoelectrons 612 — 90°. The 
phases on both panels are very similar for this orientation 
which means that the time delay of the reference photo- 
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electron does not depend significantly on its orientation 
relative to the XUV field. We have already acknowledged 
this fact when analyzing the TDSE time delay results. 

Weak sensitivity of the time delay to the field orienta- 
tion can be understood from the general parametrization 
of the DPI amplitude: 

/(fei,fe2) = [cosei + cose2]M''{Ei,E2,x) (12) 
+[cos6ii -cos6'2]A1"(£;i,^2,a;) . 

Here x = cos ^12 = cos(^2 — ^1) and the complex gerade 
and ungerade amplitudes possess the exchange 
symmetry M^^''{Ei,E2) = ±M^^''{E2, Ei). Even for 
the most severe energy sharing E2 ^ Ei, the gerade 
amplitude is still strongly dominant \A4^\ ^ [15]. 
Therefore, unless the photoelectrons are anti-parallel and 
the kinematic factor accompanying the gerade amplitude 
tends to zero, its contribution is dominant and it is 
that determine the overall DPI phase. This means that 
the timing measurement at perpendicular photoelectron 
orientation can only deliver the phase. An analogous 
measurement for would require the anti-parallel ori- 
entation which is not practicable for individual photo- 
electron streaking. However, since the relative phase of 
M.^/" can be determined independently, knowing the Ai^ 
phase will immediately deliver the missing phase of A^". 

On the top panel of Fig. 3, we show a cumulative 
phase plot argf{Ei,E2) for various combinations of the; 
reference and spectator electron energies taken at the 
perpendicular orientation 612 = 90°. The raw CCC 
data, marked by the points, arc only available for the 
E2 < Ex. To obtain the phases across the whole excess 
energy range, we fit the raw CCC data with a rational 
function and continue it analytically past the mid excess 
energy point. The energy derivative of this function, cal- 
ibrated in units of time delay, is presented on the bottom 
panel of Fig. 3. The values of time delay obtained from 
the solution of TDSE are marked by the points colored 
accordingly to match the combination of Ei and E2 from 
the CCC calculation. We observe that for these partic- 
ular combinations of the reference and spectator photo- 
electron energies, the TDSE and CCC time delays are 



quite close. 

The time delay of the reference photoelectron varies 
very rapidly with its energy but depends much weaker 
on the energy of the spectator electron. When the en- 
ergies of the both electrons are low, the time delay is 
particularly large reaching few hundred of attoseconds. 
In the opposite limit of large reference electron energy, 
the time delay becomes small. More importantly, it does 
not vary significantly with the spectator electron energy. 
Physically, this regime corresponds to the shake-off mech- 
anism of DPI in which the fast photoelectron absorbs the 
whole of the photon energy and angular momentum and 
the slow photoelectron is subsequently shaken off into 
the continuum [16]. The fast photoelectron leaves the 
atom without any significant delay, but emission of the 
slow photoelectron is delayed considerably. This delay 
becomes particularly large when the energy of the both 
photoelectrons is small and they are able to interact for 
a long time. In this regime, the main mechanism of DPI 
is the knock-out process in which the primary photoelec- 
tron impinges on the ion and knocks out the secondary 
electron into the continuum. 

In conclusion, we perform the timing analysis of the 
two-electron emission from the He atom which is sub- 
jected to a very short XUV pulse. We employ an explicit 
time-dependent treatment of the DPI process by seeking 
solution of the TDSE. We complement this procedure by 
the LOPT treatment which allows us to connect the time 
delay with the energy dependent phase of the DPI ampli- 
tude, the latter being evaluated within the CCC method. 
This opens up a possibility of a complete DPI experi- 
ment in which both the magnitudes and phases of the 
symmetrized DPI amplitudes can be determined. Such 
an experiment will require an attosecond streaking mea- 
surement on one of the two photoelectrons which can be 
performed in a close to perpendicular orientation of the 
photoelectrons. To our best knowledge, except for a very 
recent report [17], this is the first practical attosecond 
streaking measuring scheme suggested for a DPI process. 
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